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Abstract 

The magnetic deformation of the Ising Model, the thermal deformations of both the 
Tricritical Ising Model and the Tricritical Potts Model are governed by an algebraic 
structure based on the Dynkin diagram associated to the exceptional algebras En 
(respectively for n = 8,7,6). We make use of these underlying structures as well 
as of the discrete symmetries of the models to compute the matrix elements of the 
stress-energy tensor and its two-point correlation function by means of the spectral 
representation method. 



1 Introduction 



Important progress has recently been achieved in the computation of correlation functions 
for integrable models, defined either as lattice systems or as continuum theories. In addi- 
tion to the well established results on the spin-spin correlation function in the Ising model 
away from the critical temperature [|I], ^ ^ correlation functions of several important 
statistical integrable models have been obtained by means of different techniques, such 
as those discussed in the references for instance. Furthermore, for those models 

which present relativistic invariance and for which the exact S'-matrix is known, a pow- 
erful method to compute the correlation functions is provided by the Form Factor (FF) 
approach originally proposed in . This approach has proved to be extremely efficient 
because it leads to fast convergent series for the correlators, as confirmed for instance in 



121 , pr3| , |T4| . One of the most remarkable results achieved by means of the FF approach is 
the solution of the long-standing problem of the computation of the spin-spin correlator 
of the Ising Model in a Magnetic Field at T = (IMMF) |T^. The aim of this paper 
is to extend the analysis of reference to two statistical models which are very closely 
related to the IMMF, namely those relative to the thermal deformation of the Tricritical 
Ising model (TIM) and the Tricritical 3-state Potts model (TPM). The dynamics of all 
these systems are ruled by an algebraic structure related to the exceptional algebras En. 
In fact, the magnetic deformation of the Ising model highlights its underlying Eg structure 
as well as the thermal deformation of the TIM and of the TPM, which highlights respec- 
tively the Ej and the Eq structures of these models. In this paper we are concerned with 
the determination of the FFs of one of the most important fields of the above mentioned 
models, i.e. the stress-energy tensor T^y{x). The matrix elements of this operator are 
particularly simple to determine for several reasons. Firstly, all the matrix elements of 
the components of this tensor can be expressed in terms of the FFs of just one scalar 
operator, namely its trace Q{x). This operator, in turn, is related to the operator 
which deforms the conformal action by the relationship 

e(x) = 27rA (2 - 2A$) $(0;) , (1.1) 

where A$ is the conformal dimension of the field $(x). The stress-energy tensor is 
furthermore characterized by its conservation law d^T^y{x) = and by its obvious relation 
to the total energy and momentum of the system. These additional pieces of information 
are sufficient to uniquely identify the operator Q{x) and then, to reconstruct its correlation 
functions^. Some simple tests, based on certain sum rules, may be performed to check the 

^Identification of other local operators of an integrable theory may be in general more difficult. Al- 
though we are still lacking a general answer to this problem, nevertheless the series of papers [^0|, 
already yields some remarkable theoretical results. 
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efficiency of the approach. The first check involves the zeroth moment of the two-point 
correlation function G{x) = (0(x)6(O)), a quantity which is related to the amplitude 
U of the free energy per unit volume fs ~ —Urn^. Conversely, this amplitude can also 
be obtained independently by means of the Thermodynamics Bethe Ansatz ||19| , pO[] . The 
second check is based on the second moment of G{x) which is proportional to the difference 
of the central charges Cuv and Cj^ of the conformal field theories arising in the ultraviolet 



and infrared limits |21|. An important point illustrated by these comparisons is that the 
spectral representation series of the correlation function is saturated with a high degree 
of accuracy by its first terms. This implies that the analytic efforts needed to obtain an 
accurate determination of the correlation function are simplified enormouslyfj. 

The paper is organized as follows: in the next section we briefly review the bootstrap 
approach to integrable models based on the exact ^-matrix formulation. In Section 
3 and 4 we compute the first Form Factors of the stress-energy tensor of the perturbed 
Tricritical Ising Model and Tricritical Potts models respectively (the computations relative 



to the Ising model in a magnetic field may be found in the original reference [13|). Our 



conclusions are presented in Section 5. The paper also contains two appendices: in the 
first one we gather useful mathematical formulas used in our computations, while the 
second one contains the three-particle FF results of the TPM. 



2 General Features of Integrable Quantum Field 
Theories 

In order to clarify the basic principles of the bootstrap approach to the solution of an 
integrable Quantum Field Theory and to set up the notation, in this section we endeavor 
to present a concise picture of the conceptual framework on which such approach is based 
and its most important consequences. Our starting point is the theory of integrable 



scattering processes, as originally developed in [22 



2.1 Perturbed CFT and Factorized Scattering Theory 

The statistical systems analized in this article are particular integrable massive deforma- 
tions of the first minimal unitary models of Conformal Field Theories (CFT), namely the 
Ising model, the Tricritical Ising Model and the Tricritical 3-state Potts Model. The Ising 
Model is deformed along the magnetic direction whereas the deformation of the other two 
models is along their thermal direction. The three models belong to the minimal unitary 

^It is worth to notice that this remarkable behaviour of the spectral series has been also observed for 
massless models 
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series of CFT and therefore, at their critical point they may be simply realized in terms of 
a coset contraction based on the affine SU{2) algebra |2^. However, it is well known that 



at criticality they also admit an alternative realization in terms of the coset constructions 
{En)i (8> {En)i I {En)2 bascd on the exceptional algebras En- the Ising Model is associated 
to the exceptional algebra Eg, the Tri critical Ising Model to E-i and the Tricritical 3-state 



Potts Model to Eq. As shown in [g3| |29|, ^ the advantage of considering such 
alternative realizations of the three models at criticality becomes clear once they are per- 
turbed away from the fixed points along particular directions in the scaling region. Since 
we will have the opportunity to come back and comment on this point in the next sections, 
where we consider in more detail each model, here we prefer to keep our discussion as 
general as possible and focus our attention on the common features of a typical integrable 
perturbed CFT. For such a theory, the action can be formally written as 

A = AcFT + g J d^x^{x) , (2.1) 

where ^{x) is the relevant operator that moves the system away from criticality. It breaks 
the original conformal invariance and therefore introduces a mass scale in the theory. For 
integrable deformations, the dynamics of the system is supported by the presence of an 
infinite number of conserved charges (this may be argued by means of the counting argu- 
ment proposed in One of the consequences of the quantum integrability is that the 
scattering processes of the massive excitations of the theory are completely elastic. Since 
production processes are absent in the dynamics of such integrable QFT, their general 
S'-matrix element is diagonal in the number of external [in and out) particles involved 
in the collision and moreover is completely factorized into the product of two-particle 
amplitudes. Hence, to know the on-mass-shell properties of such class of QFT we have 
simply to determine the two-particle S'-matrices. For the class of models which we will 
consider in this paper, there is an additional simplification, i.e. all particles can be unam- 
bigously distinguished on the basis of their different quantum numbers (this is certainly 
the case for theories with a non degenerate mass spectrum). Under this circumstance, the 
scattering processes are completely diagonal and the two-particle scattering amplitudes 
are simply defined by the equation 

I AaiOa) A,{e,) )out = SabiOab) \ A„(^^„)^(0,) )i, , (2.2) 

where 6ab = Oa — Ob and we have used, as usual, 9a to parametrize the dispersion relation 
of the particle Aa, i.e. 

{E,p) = {maCosh6a,maSmh.6a) . 

In terms of 6ab, the Mandelstam variable s reads s = m\ + ml + 2 mam;, cosh . As 
functions of the variable 9ab (from now on simply denoted by 6*), Sab{d) are analytic 
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functions, with possible poles on the imaginary axis < ^ < ivr. They satisfy the 
functional equations 

Sab{9) Sabi-9) = I , 

(2.3) 

expressing the unitarity condition and the crossing symmetry of the theory, respectively. 
The general solution of (|2.3| ) can be written in terms of an arbitrary product of the 
functions 

_ sinh|(g + z7ra) 
smh 2(^7 — iTxa) 

where — 1 < a < 1. The parameter a is related to the position of the pole of Sa{0) which 
is located at 6' = ma. For those models with a non-degenerate mass spectrum (in which 
all particles are self-conjugate), the functional space of the solution of (|2.3| ) is instead 
spanned by an arbitrary product of the crossing symmetric functions 

. X / X tanh ^(0 + ZTTO;) 

UO) ^ s^{e)sS7, -6)= i) . . 2.5 

tanh — tTcaj 

The simple poles of are located at the crossing symmetric points 6 = ma and 9 = 
in{l — a). 

For a theory with a degenerate mass spectrum the general expression of the two- 
particle S-matrix element may be written as 

Sab{e) = n saiey- , (2.6) 

while for a theory with a non-degenerate mass spectrum we have 

Sabie) = n U^^Y" ■ (2-7) 

aeAab 

In both cases, the exponents pa denote the multiphcities of the corresponding poles iden- 
tified by the indices a. For those models which present an underlying algebraic structure 
related to a Dynkin diagram, as for instance the three models investigated in this paper, 
the labels a G Aab are integer multiples of 1/h where h is the Coxeter number of the 
associated Lie algebra. The complete two-particle S'-matrices of the Tricritical Ising and 
Potts Models in a thermal perturbation (originally computed in ^) are reported 



for convenience in Tables (2) and (7) (the one relative to the IMMF may be found in the 
original reference ||23|| ) . 

The two-particle elastic S'-matrices considered in this paper present quite a rich pat- 
tern of poles in the complex 9 plane. According to the analysis carried out by several 
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groups [0, |3^, |39|, the simple and higher odd-order poles are associated to the pres- 
ence of bound states appearing as intermediate virtual particles in the scattering processes. 
Viceversa, all even-order poles are simply due to intermediate multi-scattering processes 
with no one-particle singularities. For a simple pole of Sab{0) at 9 = iu'^f^ corresponding 
to a bound state A^, we can compute both the mass of the bound state by means of the 
equation 

ml = ml + ml + 2mamb cosm^^ , (2.8) 

and the on-mass-shell three-point coupling constant FJ^^ by taking the residue on the 
pole 

-z \imje-zu:,)s,b{e) = {T:,y . (2.9) 

ao 

In the case of a double pole of the S'-matrix placed at 6' = in a, the associated residue is 
given by 

- ^ hm {6,, - i^fS^m = (F:,F^,-)2^,e(^7) , (2.10) 

where 7 = vr — u'^ — Mgj. It is beyond our scope to discuss here in more detail the multi- 
scattering interpretation relative to the analytic structure of the S'-matrix. The interested 
reader may consult the original literature quoted above for a complete account of the pole 
structure of the S'-matrix. Some examples will be however provided in the next sections 
to enlighten the analytic structure of matrix elements of local operators. 

2.2 Correlation Functions and Form Factors 

One of the approaches that has proved to be extremely efficient in the computation of 
correlation functions for statistical models away from criticality is the spectral represen- 
tation method. For integrable models, this approach has been originally proposed in [0, || 
and further analysed and applied by different groups [7-18]. The simplest example which 
illustrates this approach is given by the computation of the two-point functions (higher- 
point functions being determined similarly). In the spectral representation approach, the 
two-point function of a local operator ^{x) may be expressed as 



C Mr? ri r? 2 

(^(^MO)) = E / , , -^n^\FL..aS(^i^---M e-l^'lS-i-'=^-^^^- , (2.11) 

n=0 ^••■>"n l^''/ 

where 

i^r„...,a„(^l,---,^n) = (0|^(0) K(^i)---<(^n)) • (2.12) 

The above matrix elements are the so-called Form Factors (FF) and, as we will briefly 
discuss below, their computation can be performed once the exact S'-matrix and the 
bound state structure of the theory are known. Inserted into (|2.11|), they give rise to 
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fast convergent series both in the infrared region ( |x| — * oo), with a corresponding 
exponential decay, and in the ultraviolet limit (|x| —* 0) where the correlation functions 
show power-law behaviours. 

The FF satisfy the so-called Watson equations, given by 

Fai,...,ai,ai+r,...,aA^l^ ■ ■ ■,0i,9i+i, . . . , 9n) = 

= Sa„a,+ A9i - 9i+l) ^a^i,...,a,+i,a.,...,a„(^l' • • • ' ^i+l, 9i, ■ ■ ■ , 9n) , 

K,a,,...,a. (^1 + 2vr^, 9,,..., 9^) = {92, • • • , ^^i) . (2.13) 

Among the solutions of these equations, there are those (called minimal solutions) charac- 
terized by the property that they have neither poles nor zeros in the strips Im9ij G (0, 2tt). 
By using the factorization properties of the underlying scattering theory, the minimal 
solution associated to a generic FF may be easily expressed in terms of the minimal 
two-particle FFs F™"(^) by 

K^Z,...,aS(^U 92,. ■ .,9^)= n Ka-i(^^ - • (2-14) 

l<i<j<n 

where the explicit expressions of F^^^{9) are given for theories with degenerate mass 
spectrum by 

FT{9) = (-zsinh|) n ^^{97" ' (2-15) 
while for theories with non-degenerate mass spectrum by 



'(^^)= (-^sinh^) n 9a{9Y'' ■ (2-16) 



The definition and the properties of the functions ha{9) and ga{9) are collected in Ap- 
pendix A. 

The minimal expression of the FFs does not carry any dependence on the specific 
operator we are considering, as it must be, since the monodromy properties derive from 
the S-matrix alone. To characterize the different operators and to take into account the 
dynamical pattern of bound states of the theory, let us consider in more detail the analytic 
structure of the FFs, starting our discussion from the occurrence of their poles. Their 
pattern may be very complicated for the multi-scattering processes of the theory. There 
are however two classes of simple order poles in the FF which have a simple and natural 
origin [Q. The first class is that of kinematical poles relative to particle-antiparticle 
singularities at the relative rapidity 9 = in with the corresponding residue given by 

- I lim^ F-,^a,au...,aA9 + i'K,9,9^,...,9^)=l^-f{ S^a, {9 - 9,)j F,,_a^ {9,, ...,9n) . 

(2.17) 
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The second class of simple order poles of the FFs which admit a simple explanation is 
that associated to bound state singularities. Namely, whenever Aa{6a) and Af,{6i,) form a 
bound state Ac{6) for the value 6ab = iW^b of their relative rapidity, then all the matrix 
elements F^^^^^ ^^{9a,0f,,0i,...,9n) involving the two particles Aa{Oa) and Ab{9f,) will 
have as well a simple order pole at the same position, with the residue ruled by the 
on-shell three-point coupling constant V^j^, i.e. (see Figure 1) 

-Z lim^(0,,-Z<,)F,V,...,a„(^a,^^6,^^l,...,^n)= r^,F^.,,...,.J • (2.18) 

In addition to the above classes of simple poles, the FFs may present poles of higher order 



relative to the underlying multi-scattering processes, as recently clarified in ||T3[| . 

A key point to understand the rich analytic structure of the matrix elements is to 
initially^ analyse the two-particle FFs. Following the analysis of ||13| (see also |]16|)) the 
two-particle FFs can be conveniently written as 

= QM^J^ , (2.19) 

where Dab takes into account its poles structure and is a polynomial in cosh 9 which 
carries the dependence on the operator yj. 

The polynomials Dab{9) are determined from the poles of the S'-matrix. The analysis of 
ref. fl^ gives the following simple rules for determining them in the case of non-degenerate 



theories: 



Dab{9)= n (^a(^))"(Pl-a(^)) , (2.20) 

aeAab 



= n + 1 , ja = n , if Pa = 2n 
ia = n , ja = n , a pa = 2n 

where Aab and pa are defined in eq. (|2.7|) . The functions 

cos 7Ta — cosh 9 
2 cos^ 



(2.21) 



= ^V7Z-r^ (2-22) 



2 



give a suitable parametrization of the pole at ^ = ivra. The above prescription can be 
also generalized to degenerate theories. In fact, referring to equation (|2.4|), one can write 



^ab{9) = n [Pc.{e)) " , (2.23) 

aeAab 



The reason is that, by factorization, FFs with higher number of particles inherit their pole struc- 
ture from the analytic structure of the two-particle channels. Moreover, the two-particle FFs play an 
important role in the theory since they provide the "initial conditions" needed for solving the recursive 



functional equations ( |2.17| ) and (|2.18| ) 
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= n + l , 


if 


Pa = 




= n , 


if 


Pa = 




= n , 


if 


Pa = 



2n + 1 s-channel pole; 
2n + 1 t-channel pole; (2.24) 
2n . 

where it is convenient to distinguish between poles associated to the direct s-channel and 
those relative to the crossed t-channel. As we will show in the sequel, the above rules 
play an essential role for implementing the boostrap program for the computation of the 
FFs in the TIM and in the TPM. Let us quote at this point the equations which will be 
often employed in the next sections. Those are: (a) the residue equations at a simple 
order pole that corresponds to a bound state 

- 2 lim^ [6 - iul,)F^,{e) = r:,F:f , (2.25) 

ab 

(see Figure 1); (b) the residue equations relative to a simple order pole induced by a 
double pole in the S'-matrix 

- z hm {9ab - iv)FaM = r^,r^-,F,e(z7) , (2.26) 

where 7 = vr — — (see Figure 2) and finally, (c) the residue equations relative to 
a double order pole induced by a third order pole in the corresponding S'-matrix (see 
Figure 3 where ip = wf^) 

hm (Oat - tui,fF,,ieat) = tK^Tl^ hm (6,, - iul)F,,{e,,) = -V'^^.Tlfi^Ff . (2.27) 

After having considered the pole structure of the two-particle FFs, let us concentrate our 
attention on the polynomial Qabi^) the numerator of ( |2.19| ). In contrast to Dab{0), 
which is only fixed by the S-matrix singularities, the polynomials Qabi^) depend, on the 
contrary, on the operator (f{x) and may be used to characterize it. An upper bound on 
the maximal degree of the polynomials Q'^bi^) has been derived in |TB]. Briefly stated. 



the argument consists in looking at the large energy limit of the FF and relating it to the 
conformal properties of the corresponding operator '^{x). Denoted by its conformal 
weight and by y^p the real quantity defined by 

hm FJ^,...,,J^i,...,0„,) -e^-l^^l 



we have 13 



< . (2.28) 

Taking into account the degree of the factor F^''"- /Dab{0) in the two-particle FF ( p.l9| ) 
by also using eq. ( |A.8| ), it is easy to translate the inequality ( p.28|) into an upper bound 
on the degree of the polynomial Qabi^)- 
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In what follows we will consider the FFs of the trace G(x) of the stress-energy tensor. 
In this case we have additional constraints for the corresponding polynomial Qab{(^)- In 
fact, the conservation law d^T^j_y{x) satisfied by the stress-energy tensor implies that the 
FFs of the trace Q{x) must contain the kinematical polynomial = {pi + ■ ■ ■ + PnY, 
with the exception of the FFs with two identical particles. For the two-particles FFs, this 
property can be expressed by means of the following factorization 

Ql{e) = cosh^ + ' Patie) , (2.29) 



2mamt 



where 



Pabie) = Y.(^tbCosh'9 , (2.30) 

k=0 

The degree Nab in (|2.30|) may be determined by implementing the inequality (|2 . 28| ) . In this 
way, the problem is reduced to determine the coefficients a^j, of the polynomials Pab- This 
goal can be achieved by applying the residue equations together with the normalization 
condition of the two-particle FF, which are expressed by 

Fgitrr) = 27rml (2.31) 

The above conditions prove in general sufficient or even redundant in number, to fix all 
the coefficients of the polynomial Pab{0). 

As mentioned in the introduction, there is now a strong evidence that the spectral 
series based on the FFs are fastly convergent. For the correlation functions of the stress- 
energy tensor, one way to test this convergence is to employ two sum rules satisfied by 
the moments Aip of the two-point function of 6(a;), defined by 

Mp = j d^x\x\P {Q{x)Q{<d)) . (2.32) 

The ffist sum rule is relative to the bulk free energy / ~ —Um^, where the amplitude U 
is related to the zero-moment A^o by 

U=T^^^o, (2.33) 
with m the lightest mass of the theory. The second sum rule relates the second moment 



to the central charge c of the original CFT, according to the formula |21 

c = . (2.34) 

47r 

By inserting the spectral representation of the two-point function (0(a;)0(O)), both 
moments can be expressed as a series on the number of the intermediate particles with 
an increasing value of their center-of-mass energy. In this way we can test the fast 
convergence of the spectral representation by comparing the truncated values of the series 
of the moments Aio and A^2, with the known value of the central charge c and the exact 
value of U computed by means of the Thermodinamic Bethe Ansatz [ll9| , pO |. 
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3 Form Factors of the Energy Operator for the 
Thermal Perturbation of the Tricritical 
Ising Model 

The Tricritical Ising model is the second model in the minimal unitary conformal series 
with central charge c = 7/10 and four relevant fields [EB]. The microscopic formulation 
of the model, its conformal properties and its scaling region nearby the critical point 



have been discussed in several papers (see, for instance ||30|, 0, |32|, |3^, ^). In the 
following we give a short review of the features of the TIM which are most relevant to the 
FF approach to integrable massive models. 

3.1 Generalities of the TIM 

The Tricritial Ising model may be regarded as the universality class of the Landau- 
Ginzburg $^-theory 

L = (V<f )2 + (76$' + ^74$' + ^73$' + 92^^ + 91^ (3.1) 



at its critical point 9i = 92 = 93 = 94 = [37]. This Lagrangian describes the continuum 



limit of microscopic models with a tricritical point, among them the Ising model with 



annealed vacancies, with an Hamiltonian given by ^ 



H = -f5 (^i'^jtitj - fi^U . (3.2) 

<ij> i 

(3 is the inverse temperature, /i the chemical potential, ai = ±1 the Ising spins and = 0, 1 
is the vacancy variable. The model has a tricritical point (/9o,/io) related to the sponta- 
neous symmetry breaking of the Z2 symmetry. At the critical point (/5o, fio), the TIM can 
be described by the following scaling fields: the energy density e{z,z) with anomalous 
dimensions (A, A) = (^, ^), the vacancy operator or subleading energy operator t{z,z) 
with (A, A) = (|, |), the irrelevant field e" with (A, A) = (|, |), the magnetization field 
(or order-parameter) a{z,z) with (A, A) = (^, J,), and the so-called subleading mag- 
netization operator a{z,z) with anomalous dimensions (^, ^). With rescpect to the Z2 
symmetry of the spin model, the spin operators are odd while the energy operator, the 
vacancy operator and the irrelevant field e" are even. 

A peculiar feature of the TIM is the presence of another infinite dimensional symmetry 
in addition to the Virasoro algebra, i.e. a hidden VT-algebra based on the Ej root system. 
This is related to the equivalent construction of the TIM in terms of the coset model 



(-£^7)1 ® {E'i)i/{E'i)2. Let us briefiy recall the coset formulation at the critical point p7 |. 
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From the theory of Kac-Moody algebras, the central charge of a CFT constructed on an 
affine Lie algebra G at level k is given by 

k\G\ 

cg = , , , , 3.3 

k + ho 

where | G | is the dimension of the algebra and he the dual Coxeter number. The unitarity 
condition for the CFT restricts the highest weight representations | A ) which can appear 
at the level k. Denoting with u the highest root, the allowed representations | A ) at the 
level k must satisfy 

'-^<k, (3,4) 

and their dimension is given by 

k + he 

where C\ is the quadratic Casimir in the representation {A}. Using a subgroup H G G, 
one can construct a CFT on the coset group G/H, with a central charge equal to 

cq H = cg-ch = ^^^^ ^ - ^^1^1 (3 6) 

Its representations are simply obtained by the decomposition of the Hilbert space 

\cgAg)= ®k [ \cg/h.^g/h)®\chAh) ] ■ (3.7) 

In the case of the TIM, /i = 18 and eq. ( p.6| ) gives c = ^. At level /c = 1, the possible 
representations are the identity 1 and the representation Ilg with scaling dimension and 
I respectively 

(EOi ^ {l,n6} = {0,|} . (3.8) 
Their components (ni, n2, ■ ■ ■ , n-j) {rii integer) with respect to the simple roots of Ej are 



given by 

1 ^ (0,0,0,0,0,0,0) 
He ^ (0,0,0,0,0,1,0) 

At the level k = 2, the representations are given by 



(3.9) 



(Erh - {i,ni,n2,n5,n6} = {0,^,§,|,|} , (3.10) 

with the corresponding fundamental weights 

Hi ^ (1,0,0,0,0,0,0) , 

U2 (0,1,0,0,0,0,0) , (3.11) 

^ (0,0,0,0,1,0,0) . 
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Hi is the adjoint representation. Using eq. ( p.7|) , the scahng dimensions of the TIM are 
recovered by the decomposition 

(o)i X (o)i = [(o)t/a/ ® (0)2] + [{johiM ® (ni)2] + [{^)tIM ® (n5)2] , 

(0)l X = [(^)t/M ® (n2)2] + [ii^hiM ® (n6)2] , (3.12) 

(!)lX(|)l = (iW®(0)2 . 

The off-critical perturbation considered in this paper is the one given by the leading 
energy operator e{z,z) of conformal weights (^j^, j^^. Note that this operator is associated 
to the adjoint of Ej. According to the analysis of this leads to a structure of the 
off-critical system deeply related to the root system of E7, as we briefly recall in the 
following. 

First of all, the off-critical massive model shares the same grading of conserved currents 
as the Affine Toda Field Theory constructed on the root system of E7, i.e. the spins of 
the higher conserved currents are equal to the exponents of the Ei algebra modulo its 
Coxeter number h = 18, i.e. 

5 = 1,5,7,9,11,13,17 (mod 18) . (3.13) 

The presence of these higher conserved currents implies the elasticity of the scattering 
processes of the massive excitations. To compute the mass spectrum and the scattering 
amplitudes, it is important to observe that, according to the sign of the coupling constant 
g in ( p.l| ), this perturbation drives the system either in its high-temperature phase or in 
its low-temperature phase. While in the latter phase we have a spontaneously symmetry 
breaking of the Z2 symmetry of the underlying microscopic spin system, in the former 
phase the Z2 symmetry is a good quantum number and therefore can be used to label the 
states. In the low-temperature phase, the massive excitations are given by kink states 
and bound state thereof, in the high-temperature phase we have instead ordinary particle 
excitations. The two phases are related by a duality transformation and therefore we can 
restrict our attention to only one of them, which we choose to be the high-temperature 
phase. In this phase, the massive excitations are given by seven self-conjugated particles 
Ai, . . . ,Ar with mass 

mi = M{g), 

57r 

m2 = 2mi cos — = (1.28557..) mi , 
18 

TT 

777,3 = 2 772,1 COS — = (1.87938..) mi , 
9 

TT 

m4 = 2 mi cos — = (1.96961..) mi , 
18 

ms = 2 m2 cos ^ = (2.53208..) mi, (3.14) 
18 
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2-K , 

me = 2 ma cos — = (2.87938..) mi , 

IT 

ruj = 4 m3 cos — = (3.70166..) mi . 
18 

The dependence of the mass scale M on the couphng constant g has been computed in 

m 

M{g) = Cgi, (3.15) 



where 

C = [4 vr^ 7(1) 7(1) 7^(^)1 ^ r(^i) ^j) rf f ) = ^.745372836 . . . , (3.16) 

where 7(x) = ^ii^x) ■ '^^^ mass ratios are proportional to the components of the Perron- 
Frobenius eigenvector of the Cartan matrix of the exceptional algebra E-j and therefore 
the particles Ai may be put in correspondence with the following representations of E-^ 
(here identified by their dimensions) 



Al - 


56 , 


A2 - 


133 , 


A, - 


912 , 


^4 - 


1539 , 


A, - 


8645 , 


A, - 


27664 


A, - 


365750 



(3.17) 



The exact S'-matrix of the model is given by the minimal S-matrix of the Affine Toda 
Field Theory based on the root system of E^. It has been calculated in |31| and is 



listed for convenience in Table 2. The structure of the bound states may be written in a 
concise way by grouping the particle states into two triplets and one singlet states 

{Ql,Q2,Q3) = (^6,^3,^0, 

{K^,K2,K3) = (^2,^4,^7), (3.18) 

(N) ^ (A,) . 

The first triplet consists of the Z2 odd particles whereas the other triplet and the singlet 
are made of Z2 even particles. The "bootstrap fusions" involving [A^] and [A^, Ki] form 
closed subsets 

N-N = N, N-Ka = K, + K2 + K,, 

Ka-Ka+i = Ka + N, Ka-Ka = Ka + Ka+1 + N. 
Including the first triplet, we obtain the following algebra 

Ka ■ Qa = Qa+1 , Ka ■ Qa+1 = Qi + Q2 + Q3 , 

Ka ■ Qa-i = Qa-i + Qa+1 , Qa - Qa = A'^-i + Ka+i , (3.20) 

Qa-Qa+1 = Ka + Ka^i + N , N-Qa = Qa-1 + Qa+i. 
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It has been observed that these bootstrap fusions are a subset of the tensor product 
decomposition of the associate representations of Ej p9|] . 



3.2 Form Factors of the TIM 

After the discussion on the general features of the model, let us consider now the problem 
of computing the FFs of the operator e(x) or, equivalently, of the trace G(x) of the 
stress-energy tensor. To this aim, the Z2 parity of the model is extremely helpful. In 
fact, because of the even parity of the energy operator, we can immediatly conclude that 
its FF with a Z2-odd (multi-particle) state must vanish. In particular, the one-particle 
FFs of for the odd particles are all zero. 

To start with the bootstrap procedure, let us consider the two-particle FF relative to 
the fundamental excitation Ai 

F,^,{e) = QfAO) , (3.21) 

where 

F^'^ie) = -I sinh(^/2) g,„{d) gy,{e) , (3.22) 

and 

Dii{e) = v,„{e)Vi„{e). (3.23) 

By using the bound ( [2.2^ ), we see that the polynomial (6*) reduces just to a constant, 
which can be easily determined by means of the normalization condition ( p.31| ), i.e. = 



27rm^. Thus -Fii(^) is now completely determined and its expression can be used to derive 
the one-particle FFs F2 and F4. Indeed, the particles A2 and A4 appear as bound state 
of the particle Ai with itself, the coupling and F^^ being easily determined by the 
residue equation (|2.9| ). By using then the equation for the bound state poles of the Form 
Factors ( |2.25| ), one gets the desired result (see Table 4). 

To proceed further, it is convenient to list the Z2 even states (the only ones giving 
non-vanishing FFs of the stress-energy tensor) in order of increasing energy, as in Table 
3. After computing , F® and F®, which are obtained by means of the same technique, 
(i.e. fixing the unknown coefficients of FFs by using the simple pole residue equations), 
a more interesting computation is represented by the two-particle FF F2i{9). The corre- 
sponding S-matrix element displays a double pole and therefore, according to eq. ( ^.2U| ), 
we have 

F^M = QZie) , (3.24) 

where 

KriO) = 97/9{0) 94/9{0) 9l/,{e) , (3.25) 
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and 

D24{e) = Vj/g{e) v,„{e) Vysie) V2/m ■ (3.26) 

Taking into account the asymptotic behaviour of the FF and eqs. ( p.29| ) and ( p.30| ), we 
conclude that in this case the polynomial P24 has degree A^24 = 1 and therefore Q2i{Q) 
reads 

cosh e + ^ ' (a° + a\, cosh 9) . (3.27) 

277121714 J 

To determine the constants 0^4 and al^, we need at least two linearly independent equa- 



tions, which are provided by eq. (|2.25|) on the fusions 



(^2,^44)^^2 and (A2,A4)-^A5. (3.28) 

Both F2 and F5 are known, of course, from previous computations. In this case, the 
double pole in the S*- matrix provides a non-trivial check for the computation. In fact, we 
have the process drawn in Figure 2, with the identification 



and respectively 



a = 2, 6 = 4, d = e = 1 , 



c = 1, (f = 27r/3, 7 = 7r/3 , 



or 



c = 3, (f = n/3, 7 = 7r/9 . 

These processes give rise to the corresponding residue equations 



-z lini (^^-z27r/3)F«4(^^) =r^ir^iF«(m/3) 



(3.29) 



hm (^ - zn/S)F^M = Hirii F^,{zn/9) 



which are indeed fulfilled. This example clearly shows the over-determined nature of the 
bootstrap equations and their internal consistency. 

The next FF in order of increasing value of the energy of the asymptotic state is given 
by the lightest Z2 even three-particle state | A1A1A2 ). The FF may be parametrized in 
the following way 

Tpmin(n \ Tpmin(n \ pminfn \ /O® 
^n2[^a,t^l>,t^c)- D,,{e,,)D,2{eac)D,2{e,c) cosh ^^,, + cosh ' ^^'^^^ 

where and Z)™" are given by equations ( p.22| ) and (|3.23|) , while 



i^lT"(^)=^7l3/18(^)^77/18(^) , (3.31) 



15 



and 

Di2ie)=Vi3/is{0)Vr/is{e) . (3.32) 
We have introduced into ( p.30| ) the term 

1 

cosh 9ac + cosh 9bc ' 

to take into account the kinematical pole of this FF at 6a = 9b + m. The polynomial Q112 
in the numerator can be further decomposed as 

gf,2(^„^,,^,) = P2pe^ , (3.33) 
where is the kinematical polynomial expressed by 

= 2m\ + m\ + 2m\ cosh^af, + 2mim2(cosh^ac + cosh^^c) • (3.34) 

The degree of P®2 can be computed by means of the asymptotic behaviour in the three 
variables 6a,b,c separately. This gives the following results for Q ~ exp [xiOi] : 

Xa = Xh = 1 and Xc = 2 . (3.35) 

Hence, a useful parametrization of the polynomial P112 is given by 

Pn2(^a5 Ob, Oc) = Pq+Pi cosh Oab + P2(cOsh Oac + COsh Obc) + P3 cosh 9ac cosh 9bc , (3.36) 

where four unknown constants have to be determined through the poles of Ff^2- By using 
the kinematical pole at Oab = "^^r and the bound state poles at Oab = "^^5 "^f and Oac = 
ij^, one obtains a redundant but nevertheless consistent system of five equations in 
the four unknown pi whose solution is given by 

= -pi = ^ = -39.74991118... , p2 = -198.2424080... (3.37) 

The other FFs which we have computed correspond to the states listed in Table 3. The 
values of the one-particle FFs are collected in Table 4, while the results concerning the 
two-particle computations are encoded in Table 5 via the coefficients a^^ of the polyno- 
mials Pab{0). 

3.3 Recursive Equations of Form Factors in the TIM 

For sake of completeness, we now illustrate an efficient technique to compute multiparticle 
FFs. This is based on recursive identities which relate FFs of the type with differ- 

ent (even) numbers of fundamental particles. Once these FFs are known, those relative 
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to Z2 even multi-particle state involving heavier particles may be obtained through boot- 
strap procedure. In general this way of proceeding is the simplest one as far as FFs with 
three or more particles are concerned. In order to write down these recursive equations, 
we can adopt the following parameterization for the 2n-particles FF -Fi.i^...^i: 



I? (a a \ — T (a a \ -f^2n.<52n(a^l, • • • , a;2n) -rr -F™"('7iA:j 

(3.38) 

Here and in the following crfe(xi, . . . , X2n) represents the symmetric polynomials of degree 
k in the variables Xi = e^' defined through their generating function 



m 



Y[{x + Xk) = Y.^"'~'^A^^^---^^m) . (3.39) 

k=l j=0 

Fii and Dn are defined by ( p. 221) and ( 3.23 ) while Hn is an overall multiplicative constant 
and Qn is a symmetric polynomial in its variables. The factors {xi + Xk)~^ give a suitable 
parametrization of the kinematical poles, while the dynamical poles are taken into account 
by the functions -Dn's. 

The polynomial Q2n in the numerator can be factorized as 

Q2n{Xi, . . . ,X2n) = CriCr2n-lP2n{Xi, . . . ,X2n) , (3.40) 

since the FF will be proportional to the kinematical term relative to the total momen- 
tum which can be conveniently written as 

F = . l<J-41j 

The Lorentz invariance of the FF requires P2n to be an homogeneous polynomial with 
respect to all the Xj's of total degree 

degP2n = 4n2-5n , (3.42) 

while the condition ( p.28| ), knowing that = 1/10, imposes an upper bound to the 
degree in a single Xj, given by 

deg,^P2n<4n- 22/5 . (3.43) 

Writing down the most general expression of P2n as a symmetrical polynomial in the 
basis of the a^'s and taking into account the above conditions, one can determine the 
relative coefficients by means of the recursive equations. A first set of recursive relations 
is obtained by plugging the parametrization of J^2n into the equation of kinematical poles 
(|2.17| ); the polynomial Qn are then solution of the recursive equation 

Q2n+2{-X, X, Xi, . . . , X2n) = -i Q2n{xi, X2n) U2n{x\Xi), (3.44) 
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where the polynomial U2n is given by 

n 

U2n{x\xi) = n n (x + e-'™Xi)(x-e'™Xi)+ (3.45) 

i=l aeAi 1 
n 

-n n (x-e-*'^"x,)(x + e^™x,). 

i=l ai^Ai 1 

The overall constants if„ have been fixed to be 

(\ — n(n— 1) 
T-r / X COS (7ra/2) \ 
16 n 9aiO) .\ \ M , 3.46 

aelii sin(7ra) ^ 

with = 27rmf . Given (52n, eq. ( p.44| ) restricts the form of the polynomial Q2n+2, al- 
though these equations cannot determine uniquely all its coefficients. In fact, polynomials 
containing the kernel factor Y[ij=iixi + Xj) can be added to a given solution Q2n+2 with 
an arbitrary multiplicative factor, without affecting the validity of eq. ( p.44| ). In order to 
have a more restrictive set of equations for the coefficients of the polynomials Q2n, we em- 
ploy the recursive equations (|2.25| ). To relate T2n+2 and T2n-, we consider two successive 
fusions AxA\ — > A2 and A2AX Ai, obtaining the following equations 

Q2n+2[.-W-, X, ^X, X2,..., X2n) = (pn M (Vj X^ Q2n{x, ^s, . . . , X2n) P2n{x\Xi) (3.47) 

where 

M = 4cos(57r/18) cos(87r/18), 
cj)^= (-l)"+iexp(-i7r(10n+ 1)/18), 
(f = exp(— i47r/9), 

and 

2n 

P2n{x\xi) = l[{x - e*«'^/\i)(x - e'^"/^x,)(a; + e'^/'xi)ix + x,) . (3.48) 

i=2 

As an application of the above equations, let us consider the determination of the FF 
JF4. Taking into account eqs. ( |3.42| ) and ( |3.43| ), we can write the following general 
parametrization for P4 as 

P4(xi, . . . , 0:4) = Cl 0-^(74 + C2 o-2a4 + C3 criO'20'3 + 04(13+ C5 . (3.49) 

From (|3.44|) , knowing Q2 = cri, one gets a first set of equations on the Cj's 

C2 = 4 (2 sin(7r/9) + sin(7r/3) + 2 sin(47r/9)), 
C5 = -4 (^sin(7r/9) + sin(47r/9)), 

C4 = Cl, (3.50) 

C3 = C5 - Cl. 
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The residual freedom in the parameters reflects the presence of kernels of eq. ( p.44| ). 
Given any solution Ql, the space of solutions is spanned by 

4 

Q2 = QI + a (yi ^3 II i^i + ^j)^ "^C (3.51) 

Eq. ( p.47|) solves this ambiguity giving the last needed equation 

_ 2 4 cos(7r/18) - 11 cos(7r/6) + 12 cos(57r/18) - 8 cos(77r/18) 

~ 3 + 5 cos(57r/9) + cos(7r/3) - 3 cos(7r/9) ' ^ ' ' 

Finally one directly computes H4 from ( |3.46| ). 

The knowledge of JF4 = Fim allows us to compute through successive applications of 



( |2.25|) almost all the FFs we needed in order to reach the required precision of the FF 
expansion of the correlation function. We have used the obtained FFs to compute the 
two-point correlation function of by means of the truncated spectral representation 
(|2.11|) . A plot of ( G(x) G(0) ) as a function of |x| is drawn in Figure 5. To control the 
accuracy of this result we have tested the fast convergence of the spectral series on the 
checks relative to the first two moments of the correlation function eqs. (|2.33|) and ( p.34| ); 
the single contributions of each multiparticle state in the two series are listed in Table 
3 and the partial sum is compared to the exact known values of the central charge c 
and of the free energy amplitude U. A fast convergence behaviour of the spectral sum is 
indeed observed and therefore the leading dominant role of the first multiparticle states 
in eq. (|2.11| ) is established. 



4 Form Factors of the Energy Operator in the Ther- 
mal Deformed Tricritical Potts Model 

In this section we will consider the FF computation for the Quantum Field Theory defined 
by the leading thermal deformation of the Tricritical 3-state Potts Model (TPM). Our 
strategy will resemble the one already applied to the TIM, with suitable generalizations 
in order to deal with this theory of degenerate mass spectrum. 



4.1 Generalities of the TPM 

The 3-state Potts Model at its tricritical point may be identified with the universality 



class of a subset of the minimal conformal model Aigj [^. Its central charge is c = 6/7. 
The model is invariant under the permutation group S3. The group S3 is the semi-direct 
product of the two abelian groups Z2 and Z3, where the Z2 group may be regarded as a 
charge conjugation symmetry implemented by the generator C. For the generator Q of 
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the symmetry, we have = 1 and QC = —CQ. The irreducible representations of 
S3 could be either singlets, invariant with respect to fl {C even or C odd) or charged 
doublets. 

The off-critical model we are interested in, is obtained by perturbing the fixed point 
action by means of the leading thermal operator e{x) with conformal dimension A = 
1/7. This is a singlet field under both symmetries, C and Q. Hence, the discrete S3 
symmetry of the fixed point is still preserved away from criticality and correspondingly 
the particle states organize into singlets or doublets. The scattering amplitudes of the 
massive excitations produced by the thermal deformation of the Tricritical Potts Model 
are nothing but the minimal S'-matrix elements of the Affine Toda Field Theory based on 
the root system of Eq (they have been determined and discussed in references ^ and 
can be found in Table 7). Poles occur at values ian with a a multiple of 1/12, 12 being 
the Coxeter number of the algebra Eq. The reason of the E^ structure in the massive 
model is due both to the equivalent realization of the critical model in terms of the coset 
{Eq)i (-£^6)1/ (-^6)2 and to the fact that the leading energy operator e{x) is associated to 
the adjoint representation in the decompostion of the fields |27|. Then, once again, one 
may apply the argument of references to conclude that the massive theory inherits 
the Eq symmetry of the fixed point. 

The exact mass spectrum consists in two doublets {Ai, Aj) and {Ah, A-j^), together with 



two singlet particle states A^ and Ah pO, p6|. Their mass ratios are given by 



mi = mj = M{g) , 

rriL = 2 cos ^ = (1.41421..) m;, (4.1) 

TT 

rrih = = 2 cos — = (1.93185..) , 

TT 

vfiH = 2 COS — = (2.73205..) , 
where the mass scale depends on g as 

M{g) = Cg^^, (4.2) 

and 



C 



7 



4 vr^ 7(1) 7(^) 7(f) 7(i)J " r(l) r(|) r(|) = 3.746559718 .... (4.3) 

The above values of the masses are proportional to the components of the Perron- 
Frobenius eigenvector of the Cartan matrix of the exceptional algebra Eq and therefore 
the particles may be associated to the dots of the Dynkin diagram (see Figure 6). Hence, 
they may be put in correspondence with the following representations of Eq (identified by 
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their dimensions) 



Al - 


- 27 , 


A-i - 


- 27 , 


Al - 


- 78 , 


Ah - 


351 


Ah - 


351 


Ah - 


2925 



By introducing the alternative notation 



(4.4) 



Al ^ 


A, 


A-i^ 


Ar 


Ah 


A2 


A-h^ 


M 


Al - 


> Bi 


Ah - 





the bootstrap fusions of this model can be written in the following compact way 25 



Ai X Ai 


= A, 


+ A2 


Ai X 


= Ar 


+ 


Ai X Ai 


= B, 




Ai X Aij^i 


= B, 


+ ^2 


Ai X Bi 


= A, 


+ A2 


~Ai X Bi 


= Ar 


+ A2 


Ai X Bi^i 


= Ai+i 


B, X B, 


= B, 


+ B2 


Bi X Bi^i 


= fii 


+ B2 



(4.5) 



It is easy to check that the above fusion rules are a subset of the tensor product decom- 



positions of the above representations of Eq [39 . 



4.2 Form Factors of the TPM 

After a brief description of the model, let us turn our attention to the determination of 
the matrix elements of the leading energy operator e{x). Our strategy will be similar to 
that employed in the case of the TIM. For the TPM, however, we have a more stringent 
selection rule coming from the Z3 symmetry. Given the even parity of the operator e{x) 
and its neutrality under the symmetry, the only matrix elements which are different 
from zero are those of singlet (multiparticle) states and they are the only contributions 
which enter the spectral representation series ( p.lll) . For convenience, the first such states 
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ordered according to the increasing value of the s- variable are listed in Table 8. Because 
of the selection rules, one very soon encounters three- and four-particle states among 
the first contributions, and therefore, the computation of FFs becomes in general quite 
involved. 

Let us briefly illustrate the most interesting FF computations of this model. As far as 
one- and two-particles FFs are concerned, we just quote the result of the computations 
since they are quite straightforward and can be obtained by following the same strategy 
already adopted for the TIM; the one-particle FFs are given in Table 9, while the co- 
efficients a^j, of the polynomials Pab{,G) of eq. ( p. 301 ) are listed in Table 10. The need to 
compute several three-particle FFs suggests however to adopt a more systematic tech- 
nique based on the recursive structure of the FF. The lowest neutral mass state is given 
in this model by a doublet of conjugated particles I and I. Hence, in order to build useful 
"fundamental" singlet multiparticle FFs we have to consider recursive equations relating 
FFs of the kind ^^^{iT.) = ^fii'i ui with an arbitrary number of particle-antiparticle pairs. 
From the knowledge of F^j^j obtained as solutions of the recursive equations, we can next 
derive (by bootstrap fusion) all the three-particle FFs we need in our determination of 
the correlation function. To write these recursive equations, let us parametrize the FFs 
as 

^n{n)i(^i^ Pi, ■ ■ ■ , Pn, Pn) = ^ (4-6) 




where 

F^"(/3,-/3J= _ (4.7) 

[ F-^'^'iP^- Pr) Otherwise . 

In these expressions Xi = e^* and am is the symmetrical polynomial of degree m in the Xj's 
(the quantities Xi and Wm are analogously defined in terms of the Pi's). The two-particle 
minimal FFs are given by (see eqs. ( |2.15| ) and ( |2.23|) ) 

FiT'^iP) _ F-^'m _ -^sinh(/3/2) hy,{P) h,/,{P) hy,{P) 



Dim D-aiP) PiMP) P2/3{P) 

^riP) _ P-nP) _ h,/M hysiP) h„m 



(4.8) 



(4.9) 



Di-m D-M PidP) 

In ( |4.6| ), Hn is just a multiplicative overall factor and Q„ is a polynomial in its arguments. 
The latter is the only unknown quantity which can be computed through the recursive 
equations. The function Q„ must be a symmetrical polynomial both in the Xj's and in 
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the Xj's separately. Furthermore, it must be symmetric under charge conjugation, i.e. 
under the simultaneous exchange Xj Xj (Vi = 1 . . . n) . Hence, it can be parametrized 
in terms of products of cr's and ct's with suitable coefficients in order to guarantee the 
self-conjugacy. The factor for this set of particles takes the form 

and, correspondingly Q„ will be factorized as 

Qn{Xi,Xi, . . .,Xn,Xn) = (Wn^lCTn + Cr„-ia,„) ((Ti + Wi)Pn{Xi,Xi, . . . . (4.11) 

The Lorentz invariance of the FF requires P„ to be an homogeneous poljTiomial with 
respect to all the x's and x's of total degree 

deg Pn = 3n^ -An , (4.12) 



while the condition ( p^.28| ), knowing that = 1/7, imposes the following upper bound 



for the degree in a single Xi (xj) 

deg,^ P„ < 3n - 74/21 . (4.13) 

These conditions drastically restrict the possible form of the polynomials Qn- 

Let us write down the form assumed by the kinematical recursive equations by using 
the parametrization ( [4.6|) 



Qn+l ( X,X,Xi,Xi,..., X}2, Xfi) i X Un {x\Xi^ ) Qn ix\i 3^1 , ... , 3^ru Xfi) , (4. 14) 

where (here An = {1/6,2/3, 1/2}) 

n 

Unix\xi,Xi) = II II {x-e'^'^x.^ix-e'^^'-'^^Xi)- (4.15) 

i=l aeAi I 

n n (x-e-*™x,)(x-e-^-(^-")x,). 
The overall constant is explicitly given by: 



n(n — 1) 
2 



Hn = 2Txm] I 2tan2(7r/6)tan^(57r/12) [] ^^(0) sin(7ra) . (4.16) 



However, the equations ([4.14|) are not in general sufficient to fix all the coefficients of 



Qn+i- A more stringent constraint is obtained by using twice eq. ( |2.25| ) in relation with 
the processes II ^ I and II I. The final equations take a very simple form: 

Qn+iivy,vy,vy,vy,x2,x2,...,Xn,Xn) = (4.17) 

= -(rii)^ yyWn{y,y,X2,X2,...,Xn,Xn) Qn{y,y,X2,X2,...,Xn,Xn) , 
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where r] = e*'^/^ and 



(4.18) 



[xi + Xi){xi — e 6 xi){xi — e 6 Xi)[Xi — e 2 — e 2 Xij 



■ + Xi)(a;i + a;j)(a;i - e"6"a;i)(a;i - e 6"'a;.)(a;^ - eTa;.)(a;^ - e e^'x^ 



i=2 



Let us now illustrate how this procedure works in the case of ^2(«I)- start from Fj^j, 

using eq. ( |2.31D we easily obtain Qi = 1 and Hi = 27imf. From eqs. ( |4.12D and ( [4.1.j| ), 
the general parametrization for P2 is given by 



P2{Xi,Xi,X2,X2) = Ci ((T2 + 0-2) + C2 (o-iCTsCTi + O-iO-sO-i) + 



(4.19) 



+C3 (cr^(J2 + 0-^(72) + C4 cr^cr^ + C5Cr2a2 

Equation ( |4.14| ) gives four equations for the five parameters 

C4 = -(3 + V3), 
C2-C3 = -3(2 + V3), 
ci - C2 = 3 + 2V3, 
2c2 + C5 = -18 - IOV3, 



while eq. ( |4.17| ) solve the residual freedom yielding 

9 + 5^3 



Cl 

C3 
C4 



3(5 + 3\/3) 
2 ' 
3(1 + v^) 
2 

C5 = -(3 + VS) . 



(4.20) 



(4.21) 



Once we have determined Hi and P2, we can obtain ^2(il) from eqs. ([4.6|) and ( [l.llj ). 
From this four-particles FF it is also easy to obtain the three-particles FFs Fm, F^j^, 
Fi I h applying the residue equation (|2.18| ) at the fusion angles u^j and m~ respectively. 
The explicit expressions of these three-particle FFs are given in Appendix B. 

The FFs calculated for the TPM can be used to estimate the two-point function of the 
stress-energy tensor whose plot is shown in Figure 7. The convergence of the series may 
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be checked through the sum-rule tests: the contributions of each multiparticle state are 
hsted in Table 8 where the exact and computed values of c and U are compared. A very 
fast convergence behaviour is indeed observed which supports the validity of the spectral 
approach to correlations functions in integrable massive models. 

5 Conclusions 

In this paper we have applied the Form Factor approach to estimate the correlation 
functions of the stress-energy tensor in the Tricritical Ising and 3-state Potts models. 
Both models have been perturbed away from the critical point by means of the leading 
thermal operator. In our computation, an important role has been played by the discrete 
symmetries of the two models, a Z2 symmetry for the TIM and a Ss symmetry for the 
TPM. These symmetries have in fact selected the appropriate particle states entering 
the spectral representation of the two-point correlator of the stress-energy tensor. This 
correlator has been plotted in Figure 5 for the TIM and in Figure 7 for the TPM: these 
plots are expected to be extremely precise in the large distance region {mR > 1) and 
sufficiently accurate in the crossover and ultraviolet regions {mR < 1). Obviously, a 
definite confirmation of their validity can only be obtained by a comparison with some 
experimental data or numerical simulations. 

Acknowledgments. It is a great pleasure to thank G. Delfino for useful discussions and 
helpful comments. 
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Appendix A 



In this appendix we collect some different explicit representations of the functions ga{d) 
and ha{0) together with some useful functional relations. 

Let us start by considering the non-degenerate field theories. In this case, the basic 
functions ga needed to build the minimal form factors are obtained as solution of the 
equations 

ga{0) = -fa{0)gai-e) , 

(A.l) 

gaiiTT + 9) = ga{iT^ - 0) , 

where 

tanh + iira) 



fa{0) 



tanh I — iTra) 



(A.2) 



They are called minimal solutions because they do not present neither poles nor zeros in 
the strip ImO G (0,27r). They admit several equivalent representations. The first is the 
integral representation given by 



ga = exp 



~dtoosh|(a-l/2W^^^ . - 
t cosht/2sinht ^ ' \ 



(A.3) 



where 9 — m — 6. The analytic continuation of the above expression is provided by the 
infinite product representation 



g.io) = n 

fc=0 



1 + 



e/2-K 



/ - \ 2 
1 I ( 



U+l+f 



1 + 



1 + 



e/2-K 



k+l 



(A.4) 



which explicitly shows the position of the infinite number of poles outside the strip 
ImO e (0, 27r). Another useful representation particularly suitable for deriving functional 
equations is the following: 



ga{e) = n 



r 


[l + k + 


2 + ^ 27r y 


r 


(l + A:- 




r 


(l + k- 


f + 


r 


a+k+ 


f+4) 



t=or^ (i + ^-f) (i + ^ + C 

where we have used the notation 

r(a + i9/2n) ^ = r(a + i9/2Tr) r{a - i9/2TT) . 



(A.5) 
(A.6) 
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A representation that is particularly suitable for numerical evaluations is the mixed one 



N-l 

9a{o) = n 

k=0 



X exp 



e/2Tr 
fc+l-f 



e/2TT 



e/2n 



1 + 



e/2TT 



k+l 



(A.7) 



dt cosh \Ul-2a) 



(iV + l-iVe-2*)e-2^*sin2 



t cosh ^ sinh t 



et 

2^ 



In this formula is an arbitrary integer number which may be adopted to obtain a fast 
convergence of the integral. 

Using the integral representation ( |A.3| ), it is easy to establish the asymptotic behaviour 
of Qa 

g^{9) ~ e\'\/^ for 
The function Qa is normalized according to 



^ ^ oo. 



and satisfies 



with 



5(„(i7r) = 1 , 



9o{6) = 9i{d) = -isinh- . 



The above functions satisfy the following set of functional equations 

9a{0 + iT^)9a{d) = — i ^ (sinh^ + isinvra) , 

sin na 

\ 9a+-y[^)9a--y[^) ) 

sinh I [6' — i{a — l)7r] sinh |[^^ + i{o. + l)??] 
sin^ ^ 



gi-a{Q)goc~\{Q) 



(A.9) 

(A.IO) 
(A.ll) 

(A.12) 
(A.13) 
(A.14) 



Let us turn our attention to the field theories with a degenerate mass spectrum. In 
complete analogy with the previous case, we start our analysis from the minimal solutions 
of the equations 

K{e) = ~s^{e) K{-e) 

(A.15) 

haiif^ + 6) = haiiii — 9), 

where 

sinh - ( f) 4- iTrn^ 

(A.16) 



sinh i (61 + 27ra) 

Said) - 



sinh h{9 — ina) 
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The function ha{0) is explicitly given in terms of the following equivalent representations 



ha{d) ^ exp 



.°° dt sinh [(1 - a)t] . 2 

^ / 9 

t sinh t 



sin^(^i/27r) 



K(e) = n 

fe=0 



2\ 



. (n^fr r^(fe + i + f )r(fe + 1 - f - f )r(fe + 2 - f + g) 
Ai, r2(^ + i-f)r(^ + f-S)r(^ + i + f + f) 

The mixed representation is in this case 



(A.17) 



(A.18) 



(A.19) 



N+l 

K{e) = n 

fe=0 















2 

) 





X 



(A.20) 



X exp 



2/ -(Ar + l-7Ve-*)e-^*''"^f('-")'] 



T 



sinh i 



and the asymptotic behaviour depends on the value of a 

(l-a)|9| 

hoAp) ~ e 2 for ^ ^ 00 . 
The function ha is normalized according to 

^a(«7r) = 1 

and satisfies the following functional equations: 



(A.21) 



(A.22) 



/io(^) = -isinh(^/2), 
/ii(^) = 1, 



(A.23) 



/il+a(^)=/ir-a(^), 

The basic "composition rules" for products of /Iq's are: 

hJO + ztt) hi JO) = \ cosh 

^ ^ i «v ^ cosh(^) 2 



(A.24) 
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where the polynomial V is defined in ( p. 221) of Section 2. 

Finally, since fa{d) = Sa(6')si_Q,(6'), the function can be obtained from the haS 
simply through: 

9a{e) = K{e) h.^e). (A.25) 



Appendix B 

In this appendix we briefiy report the results of the three-particle FFs relevant for our 
computation in the TPM. These FFs have been derived by applying the residue equations 
( |2.18| ) to the four-particles FF F^^^j, as explained in section ^. In writing their final form, 
we have extensively used the formulas reported in Appendix A. The two-particle minimal 
FFs -F^*"" appearing in the expressions which follow are defined by eq. (|2.15|) while the 
Dab factors parametrizing the dynamical poles are defined by eq. ( |2.23| ). 

The FF F^i is obtained from F^^^j through the residue equation at Ujj = 2i7r/3 

F^,ie„e„e,)=(l[^^fM] (^m^ + 2m^Y.^osHe,,)] . (B.l) 



i<j Dii{6ij) 



i<3 



In this expression one immediately recognizes the "minimal" part, the dynamical poles 
and the polynomial, while the only remaining polynomial in the cosh(^jj)'s allowed by 
eq. ( p.28| ) is simply a constant given by 



a'lii = -102.3375342.. 



The FF F®^, is obtained from F^^^j by using eq. (|2.18|) , with Uj^^ = in/2. Its final expression 
is given by 



Tpmin(n \ rpmintn \ rpmintn \ 



2mf + m| + 2mf cosh(6'i2) + 2mi mL(cosh.{6is) + cosh(^^23)) 



23, 



cosh(6'i3) + cosh(6'23) 
a°j^ (1 - cosh(6'i2) + 2 cosh(6'i3) cosh(6'23)) + a,^|^(cosh(6'i3) + cosh(6'; 



23, 



This expression also exhibits a kinematical pole due to the presence of a particle-antipar- 
ticle pair II. Moreover there is a nontrivial polynomial in the cosh(6'jj)'s with coefficients 
given by 



flO^^ = -70.50661963 . . . , 
a]j^ = -235.9197474... . 
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Finally, applying eq.( p.l8|) to F^^j at = in/6 one obtains 



2 mf + ml + 2mf cosh(6'i2) + 2mi nih (^cosh(6'i3) + cosh(6'23~ 



^(^uh + (^lih (cosh(6'i3) + cosh(6'23)) + a^uh cosli(6'i2) + cosh(6'i3) cosh(6'23; 
where the coefficients af;^ are given by 



a^iih = 78134.00044 



ajif, = 72661.45729. 



ifi^ = 31793.68905. 



af^f^ = 43430.98692 
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Table Captions 



Table 1 Particle spectrum, mass ratios and Z2-charges in the TIM. 

Table 2 Two-particle S'-matrix elements of the TIM; the notation (x) = fx/h{0) has 
been followed, where = 18 is the Coxeter number of E-i and the function is 
defined in eq. (|A.2| ). Superscripts label the particles occurring as bound states at 
the fusion angles = xn/h. 

Table 3 The first Z2-even multiparticle states of the TIM ordered according to the in- 
creasing value of the center-of-mass energy and their relative contributions to the 
spectral sum rules of the central charge c and the free-energy amplitude U . 

Table 4 One-particle FFs of the Z2-even particles of the TIM. 



Table 5 Coefficients which enter eq. ( |2.3CI| ) for the lightest two-particle FFs of the TIM. 

Table 6 Particle spectrum, mass ratios and Z^-chaiges in the TPM. 

Table 7 Two-particle S'-matrix elements of the TPM. In this case [x] = Sx/h{(^), where 
Sa is defined in eq. ( |A.16| ) and /i = 12 is the Coxeter number of Eq. 

Table 8 The first Za-neutral multiparticle states of the TIM ordered according to the 
increasing value of the center-of-mass energy and their relative contributions to the 
spectral sum rules of the central charge c and the free-energy amplitude U . 

Table 9 One-particle FFs of the Z2-neutral particles of the TPM. 

Table 10 Coefficients which enter eq. (|2.30| ) for the lightest two-particle FFs of the 
TPM. 
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J. ' 7 

particle 


mass /mi 


Z2 charge 


Ai 


1.00000 


-1 


A2 


1.28558 


1 


^3 


1.87939 


-1 


A, 


1.96962 


1 


A, 


2.53209 


1 


A, 


2.87939 


-1 


A, 


3.70167 


1 



Table 1 
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a b 


^ CLu 


1 1 


2 4 

- ClOl (2) 


1 2 


1 3 

fl3'l (7) 


1 3 

J- tJ 


2 4 5 

- (U) (10) (6) 


1 4 


13 6 

(17) (11) (3) (9) 




3 6 

(14) (H) (6)^ 




4 5 7 

- (^6) (^2) (A) (^o)'^ 


1 7 


6 

fl5^ (9) (F))^ (7)^ 


2 2 


2 4 5 

(12) (S) (2) 


2 3 


13 6 
(15) (11) (5) (9) 


2 4 


2 5 

(14) (H) (6)^ 
yi^) yo) y\j) 


2 5 


(17) (13) (I) (7)"^ (9) 
yii ) yio) yo) yi ) yij) 


2 6 


3 

(15) (7)^ (5)^ (9) 


2 7 


5 7 

fieUlOl^ (4:)'^ (6)"^ 


3 3 


2 7 

- (^4) (2) fSl^ (^2)'^ 


3 4 


(15) (5)^(7)^(9) 

yj.'jj y'jj yi J y^j 


3 5 


(16) (10)3 (4)^ (6)^ 


3 6 


2 5 7 

- (16) (12)3 (8)3 (4)2 


3 7 


(17) (13)3 (3)^ (7)^ (9)^ 



Table 2 (Continued) 
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a b 




4 4 


(12) (IQ)^ (7) (2)'^ 


4 5 


2 4 7 

(15) (1^)^ (7)^ (<^) 


4 6 


1 6 

(17) (11)^ (3)'^ (5)'^ (9)"^ 
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4 5 

(16) fl4)3 (f,)^ (8)"^ 
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fl2V^ r2l2 ('412 /o\4 
yiz.) yz,) y^) yo) 


Ft fi 

(J W 
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2 4 7 
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Table 2 (Continuation) 
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state 


1 2 
s/mi 


c-series 


T T 

U-series 


A2 


1.28558 


0.6450605 


0.0706975 


Aa 


1.96962 


0.0256997 


0.0066115 


Ai Ai 


> 2.00000 


0.0182735 


0.0071135 


-4, 


2.53209 


0.0032117 


0.0013783 


A2 A2 


> 2.57115 


0.0032549 


0.0025194 


Ai A3 


> 2.87939 


0.0012782 


0.0020630 


A2 A4 


> 3.25519 


0.0003010 


0.0007277 


Ai Ai A2 


> 3.28558 


0.0007139 


0.001184 


A, 


3.70167 


0.0000316 


0.0000287 


A3 A3 


> 3.75877 


0.0000700 


0.0001173 


A2 A5 


> 3.81766 


0.0000860 


0.0001581 


partidl stun 


0.6980109 


0.0911150 


exact value 


0.7000000 


0.0942097 



Table 3 
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0.9604936853 


Ff 


-0.4500141924 


Ff 


0.2641467199 


Ff 


-0.0556906385 



Table 4 



"ll 


6.283185307 


"13 


30.70767637 


"22 


15.09207695 


^22 


4.707833688 


"24 


79.32168252 




16.15028004 


"33 


295.3281130 


'^33 
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^33 


123.8295119 


"25 
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'^25 


4062.255130 


^25 


556.5589101 



Table 5 
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parncJe 


mass /mi 


Z3 charge 


Ai 


1.00000 


g27ri/3 




1.00000 




Al 


1.41421 


1 


Au 
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g-27ri/3 
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2.73205 
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Table 6 
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[101 [81^ [61^ [41 2 [21^ 

[±UJ [OJ [UJ [^J 


h h 

IV IV 


I h 

[101 [8P [6P [4P [2P 

[lUJ [OJ [UJ [^J [^J 


h h 


H 

-[1012[812 [61^ [41^ [21 

[J.UJ [OJ [UJ [Z,J 


h H 


[A] [k [7f [5]^ [3]3 [1] 


h H 


[A] [k [7]^ [5]^ [3]3 [1] 


H H 


- [10]3 [8]5 [6]« [4]5 [2]3 
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state 


1 2 

slm\ 


c-series 


M-series 


Al 


1.41421 


0.7596531 


0.0705265 


AiAj 


> 2.00000 


0.0844238 


0.0229507 


Ah 


2.73205 


0.0029236 


0.001013 


Al Al 


> 2.82813 


0.0021119 


0.0019380 


AiA^ 


> 2.93185 


0.0023884 


0.0016745 


AjAh 


> 2.93185 


0.0023884 


0.0016745 


Al Al Al 


> 3.00000 


0.0004215 


0.0004925 


Aj Aj Aj 


> 3.00000 


0.0004215 


0.0004925 


Al Aj Al 


> 3.41421 


0.00159 


0.000251 


^hAj^ 


> 3.86370 


0.0000504 


0.0001476 


At Al Ah 


> 3.93185 


0.000089 


0.0002015 


-4i Aj .4^ 


> 3.93185 


0.000089 


0.0002015 


Al Aj Al Aj 


> 4.00000 


0.0000959 


0.000381 


partial sum 


0.8569765 


0.1019449 


exact value 


0.8571429 


0.1056624 
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1.261353947 




0.292037405 



Table 9 



-'a 


6.283185307 




21.76559237 


All 


9.199221756 




25.22648264 


-In 


414.1182423 




565.6960386 




175.0269632 
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Figure Captions 



Figure 1 Diagrammatic interpretation of the process responsible for a single-pole in a 
Form Factor. 

Figure 2 Diagrammatic interpretation of the process responsible for a double-pole in a 
Form Factor. 

Figure 3 Diagrammatic interpretation of the process responsible for a triple-pole in a 
Form Factor (here = u^ab)- 

Figure 4 Dynkin diagram of E-j and assignment of the masses to the corresponding dots. 

Figure 5 Plot of the correlation function ( B(x)G(0) )/mf versus the scaling variable 
nil \x\ in the TIM. The spectral series ( ^.111 ) includes the FF contributions relative 
to the multiparticle states in Table 3. 

Figure 6 Dynkin diagram of and assignment of the masses to the corresponding dots. 

Figure 7 Plot of the correlation function ( 0(x) 0(0) )/mf versus the scaling variable 
mi \x\ in the TPM. The spectral series ( p.ll| ) includes the FF contributions relative 
to the multiparticle states in Table 8. 
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Figure 5 
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